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1 Introduction 


The purpose of this final report is to summarize our accomplishments over the period May 
93 to October 1996 for the investigation entitled “Retrieval of soil moisture and roughness 
from the polarimetric radar response" supported by the Terrestrial Ecology Program of NASA 
Headquarter. The main objective of this investigation was the characterization of soil moisture 
using imaging radars. In order to accomplish this task, a number of intermediate steps had 
to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of 
electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing 
of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly 
influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil 
moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces 
was considered and a hybrid model that relates the radar backscattering coefficient to soil 
moisture and surface roughness was developed. This model is based on extensive experimental 
measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave 
frequencies over a wide range of moisture conditions and roughness scales in conjunction with 
existing theoretical surface scattering models in limiting cases (small perturbation, physical 
optics, and geometrical optics models). Also a simple inversion algorithm capable of providing 
accurate estimates of soil moisture content and surface rms height from single-frequency multi- 
polarization radar observations was developed. The accuracy of the model and its inversion 
algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil 
surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized 
phase difference into the model. Experimental data in conjunction with numerical simulations are 
used to relate the soil moisture content and surface roughness to the phase difference statistics. 
For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random 
surfaces was developed. Finally the scattering problem of short vegetation cover above a rough 
soil surface was considered. A general scattering model for grass-blades of arbitrary cross section 
was developed and incorporated in a first order random media model. The vegetation model and 
the bare-soil model are combined and the accuracy of the combined model is evaluated against 
experimental observations from a wheat field over the entire growing season. A complete set of 
ground-truth data and polarimetric backscatter data at L-, C-, and X-band over a wide range 
of incidence angles were collected. Also an inversion algorithm for estimating soil moisture 
and surface roughness from multi-polarized multi-frequency observations of vegetation-covered 
ground is developed. 

In what follows a summary of major accomplishments with pertinent references are provided. 


2 Summary of Major Accomplishments 


The following tasks were accomplished during the past year: 
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1. Development of a polarimetric calibration technique for distributed targets: 

Quantitative analysis of the measured backscatter from rough surfaces required the use of 
a very precise radar calibration procedure. Traditionally, calibration of distributed targets 
is performed using calibration methods developed for point targets and using the illumi- 
nation integral to calculate the backscattering coefficients. Using this method, possible 
phase variations or antenna cross-talk variations (between orthogonal polarization chan- 
nels) across the beam are totally ignored, which may compromise the calibration accuracy. 
To rectify this deficiency of existing calibration techniques, a new technique was devel- 
oped with which the radar polarization distortion matrix is characterized completely by 
measuring the polarimetric response of a sphere over the entire mainlobe of the antenna, 
rather than along only the boresight direction [1]. 

2. Development of a Semi-Empirical Scattering Model: 

At microwave frequencies many natural surfaces do not fall into the validity regions of the 
theoretical models, and even when they do, the results based on these models fail to agree 
with the measurements. To circumvent this difficulty, we took a rather unique approach 
to the problem. Using our calibrated polarimetric radar system, extensive backscatter 
measurements covering a wide range of surface roughness and moisture conditions were 
conducted. The statistical parameters of the surfaces were determined by a laser pro- 
filometer and soil moisture was measured by taking samples and using dielectric probes. 
Using the polarimetric data and the expected asymptotic behavior of backscattering co- 
efficients in limiting cases, a semi-empirical model (SEM) was developed and its accuracy 
was tested [2, 3, 18]. item Development of a fully polarimetric hybrid model for bare soil 
surfaces: 

The semi-empirical model for bare soil surfaces which was developed in the first phase 
of this investigation was improved. The improved model is made fully polarimetric by 
including the phase difference statistics [19, 20]. 

3. Development of a Level-1 SAR image classifier: 

A radar image classifier was developed for the purpose of identifying statistically homoge- 
neous distributed targets in the imaged scene. This task was important it bears directly 
on to the soil-moisture estimation objective. Basically it is mandatory to identify the 
cover class of a given pixel before any algorithm is applied to it. This classifier operates 
on calibrated L- and C-band data and segments the image into four categories: (1) Bare 
surfaces, (2) Short vegetation, (3) Tall vegetation, and (4) Urban. The accuracy of the 
classifier was tested using independent images acquired in different seasons [4], 

4. Development of a dielectric probe for measuring the soil moisture: 

Traditionally coaxial probes have been used for measuring soil moisture in fields. These 
probes are very sensitive to pressure and have a very small contact area. Accurate mea- 
surement using these probes require a large number of sample measurements which is very 
time consuming. To rectify this problem we designed and completed a prototype soil probe 
based on microstrip resonators that is not only very accurate, but also very convenient to 
use. With this new probe it is also possible to predict the soil type to some extent [5]. 
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5. Development of in situ and in vivo vegetation dielectric measurement technique: 

A waveguide technique for measuring the dielectric constant of vegetation needles was 
developed. Basically the dielectric constant of the sample is calculated from the measured 
reflection coefficient using a novel inversion algorithm based on an eigen-analysis of the 
impedance matrix of the method of moments solution [6]. 

6. Development of bistatic scattering measurement facility: 

In order to characterize soil moisture and roughness parameters of a surface covered with 
randomly oriented particles (as in the vegetation case), the bistatic scattering character- 
istics of the surface are required. The specular direction is the most important bistatic 
direction in this type of problems because the mean-field (coherent scattering) is not zero. 
For this purpose a fully polarimetric bistatic scattering system was designed and tested. 
The Bistatic Facility is currently configured to perform bistatic measurements in the 8.5-10 
GHz and 34-35 GHz frequency range [7, 8]. 

7. Development of scattering model for inhomogeneous rough surfaces: 

One major problem with existing theoretical models for rough surfaces as applied to soil 
surfaces is their inability to treat inhomogeneous dielectric profile. Recently we developed 
a theoretical model that can take into account the effect of vertical inhomogeneity. For 
the first time a complete second order fully polarimetric scattering model for stratified 
rough surfaces with relatively small rms height. The model was experimentally verified 
using careful controlled experiments [9, 21]. 

8. Development of a numerical technique for rough surfaces: 

Numerical simulations are necessary to indicate the trend of the scattering behavior, par- 
ticularly where theoretical solutions and/or controlled experimental data do not exit. In 
specific three different approaches addressing different problems in numerical evaluation 
scattering behavior of rough surfaces were considered: 

• One problem in numerical computation of scattering from rough surfaces is related to 
the fact that the scattering is usually computed from finite samples of rough surfaces. 
In order to eliminate the edge efFect of finite samples, the magnitude of the incident 
wave is tapered. In these methods, the dimension of surface samples must be chosen 
large enough so that the backscatter becomes independent of footprint size [22]. To 
rectify this problem, we developed a novel numerical method that can be illuminated 
by plane waves. To eliminate the effect of edges, tapered resistive sheets are used. 
It is shown that the new method is numerically more efficient than the traditional 
methods [12]. 

• The numerical techniques available in the literature are limited to perfectly conduct- 
ing surfaces. However, the dielectric properties of soil surfaces are very much differ- 
ent from those of perfect conductors. Due to moisture variations and inhomogeneity 
among the constituent particles, the soil medium usually behaves as an inhomo- 
geneous dielectric. The numerical simulation developed here has the capability of 
handling inhomogeneous dielectric surfaces. This simulation was used to evaluate 
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the effect of moisture profile and correlation function on the backscatter data. It 
was found that both the correlation function and the dielectric inhomogeneities play 
a very important role [10, 11]. 

• One limitation of existing numerical methods is their lack of efficiency. To improve 
the computation time wavelet basis function are used to produce spars impedance 
matrices which are easily invertible [23, 24]. 

9. Retrieval of soil moisture form vegetation covered surfaces: 

Under this task some basic studies were carried out. We developed forward scattering 
models for characterizing the dependence of radar backscatter to soil moisture in the 
presence of short vegetation. The following models were developed: 

• Conducting of field measurements for surfaces covered with short vegetation: 

A complete set of ground-truth data and polarimetric backscatter data were collected. 
A wheat field was chosen as the test field and radar measurement was conducted 
over the entire growing season. The backscatter data collected on the wheat field 
was calibrated and processed completely. This data was collected in summer 1993 
during the whole growing season (mid-April to mid-September) at 6 incidence angles 
and 3 frequencies (L-, C-, and X-band). The experimental results show a significant 
backscatter sensitivity to soil moisture for all three frequencies and under different 
vegetation biomass conditions [15]. 

• Development of scattering model for grass blades: 

A analytical scattering model for grass-blades of arbitrary cross section was devel- 
oped. Closed form expressions for the polarizability tensor elements are derived. 
This model is valid when the all cross section dimensions are small compared to the 
wavelength [13]. 

• Generalization of Rayleigh-Gans model: 

A theoretical scattering model for long, thin dielectric cylinders of arbitrary cross sec- 
tion and electrical length was developed. This is an extension of our previous model 
for electrically long cylinders. Using this solution we are able to model scattering 
from grass blades of arbitrary cross section and curvature very accurately [14]. 

• Coherent Scattering Model for Short Vegetation : 

A coherent scattering model with high fidelity was developed. In this model, the 
effects of coherence among the scatterers as well as non-uniform illumination are 
taken into account. The scattering formulation for the grass components accounts 
for fine geometrical features of the grass blades such as curvature and cross section. 
[16, 17] 
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